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ABSTRACT 


Amplitude and phase relationships between sinusoidal pressure dis- 
turbances in the exit plane of a conical nozzle and mass flux response 
were experimentally determined for the subsonic flow of air. in the 
nozzle. The experimental results were compared with values predicted 
by using the modified asymptotic WKB method. | 

Experimental measurements were made with a constant temperature 
hot wire anemometer and a condenser microphone. Disturbance frequen- 
cies from 20 to 250 Hz were investigated for Mach numbers up to .34. 

The experimental values of amplitude factor agree quite well with 
analytical predictions for all Mach numbers tested. Experimental 
values for phase angle differed from those predicted but exhibited a 


imitar trend. 
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NOMENCLATURE 


Speed of sound 

speed of sound at stagnation conditions 
flow area 

flow area where M = 1 


arbitrary constants 
substantial derivative 


defined in Eq. (24) 

defined in Eq. (29) 

enthalpy | 

stagnation enthalpy 

Va 

Mach number 

defined in Eq. (62) 

amplitude factor 

Pressure 

velocity (magnitude) 

velocity vector 

radius measured from apex of cone 
radius where M = ] 

entropy 

modified Strouhal number wr/q 
temperature 

Stagnation temperature 

time 


real and imaginary components of Z, 





Oo oO 


<| 


real and imaginary components of i. 
hot wire voltage 

hot wire zero voltage | 

frequency number wh/a 

= Mc 
defined in Eq. (65) 
defined in Eq. (34) 


ratio of specific heats (= 71.4 for air) 
r/R | 

defined in Eq. (25) 

(pq) /9q 

transformed Z% 

p’/p 

density 

Mt netStonTecs time a t/R 


dimensionless enthalpy perturbation H’/a,* 


dx 
dé 
phase angle 


angular frequency 
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I. INTRODUCTION 


The engineering application of nozzles in the control of fluid flow 
is extremely widespread. It ranges in scope from small fluidic devices 
to large rocket thrusters. In the steady state there exists, for a 
given absolute pressure ratio across the nozzle, a corresponding mass 
flow rate. The introduction of pressure fluctuations at the nozzle 
boundaries affects the flow and certain dynamic characteristics are ex- 
hibited for a given set of conditions, i.e., nozzle geometry, type of 
fluid medium, and the nature of the pressure fluctuations. Because 
such fluctuations are quite common, especially in systems possessing 
cyclic operating SeponeHiSs. a Knowledge of the dynamic characteristics 
becomes important. In the past, a quasi-steady relationship between 
mass Tiow rate and pressure fluctuations was assumed. Thus, an unsteady 
Situation was reduced to a sequence of instantaneous steady state cor- 
respondences. The validity of this assumption is questionable. 

Oppenheim and Chilton [3] conducted a literature survey on pulsating 
flow measurement. Most of the investigations noted had to do with the 
effect of pressure fluctuations on the time average of flow. Elrod [1], 
working in connection with a stability analysis of externally pressurized 
thrust bearings, performed a frequency response analysis to predict the 
amplitude and phase relations between the instantaneous flow rate and a 
Sinusoidally fluctuating exit pressure in an infinite conical nozzle. 
Although the chosen geometry was extremely simple, Elrod [1] stated 
"It is believed that the essential characteristics of a compressible 
fluid expanding and accelerating from a stagnation reservoir are well 
enough realized for the results to have semi-quantitative practical 


significance." He further assumed that the flow was adiabatic, 
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frictionless, and irrotational. The analysis using the WKB method is 
asymptotically valid for high frequency number W=wR/a and low Mach 
number. Unfortunately, the numerical results were wrong because of 
a wrong choice of sign for the divergence term of the developing equa- 
tions. 

Chiang, et al., 2s realizing the generality of Elrod's analysis 
(it can be applied to either converging or diverging flow), corrected 
the choice of sign and extended the region of validity to lower fre- 
quency and higher Mach number by modifying the WKB method. They also 
investigated the problem of pulsating flow in a finite conical nozzle. 

The objective of this thesis was to experimentally investigate the 
mass flux - pressure relationship for pulsating subsonic flow ina 
conical nozzle and to compare experimental results with values pre- 


dicted using the modified WK8 solution. 


II. THEORETICAL ANALYSIS 


Elrod [1], in his original analysis, assumed an adiabatic, friction 
free (i.e., isentropic) flow originating in a stagnation reservoir. 
Under these conditions, the flow remains irrotational in accordance 


with Kelvin's theorem and the following governing equations pertain: 


Momentum 
ba+4LvyPp=o0 (1) 
Dt 
Continuity 
=O (2) 


meee’ % 





which can also be written as 


A +V-G=0 (3) 


Gibb's Equation for a simple compressible substance 


Tds-= gh-— dp = 6 (isentropic) (4) 


For irrotational flow: 


a 


BG ZV GG) = VE ° 
(see Robertson [4], p. 46.) | 


Usang Eq. (5), Eq. (1) can be rewritten as 
I: y oy ; 
o& LV P- UF (6) 
Stagnation enthalpy is defined as 


, 2 
H=A+— , (7) 
2. | 
Equations (4), (6), and (7) combine to yield 


=e 


2% =~ 7H | (8) 


r 


The scalar product of q with Eq. (1), together with Eq. (4), gives 


» 


Dig i = --oOe 2 Dae an. 9 
-" = ea arn * De at Ti 





which can also be written as 


DH. | ap 


pt ¢ ot (10) 


by considering a time derivative form of Eq. (4). 


The acoustic velocity is defined by 


(? a 


Therefore, by the chain rule, Eq. (10) becomes 


I DH. 1 20 . abne | 
Spe @ at ee (2) 


n|- 


Taking the total derivative of Eq. (12) gives 


Df. be |= 2 Ln p ae 
| | 


De Be Oa 

but 

Be Bt RTO EVE 
de (ort) “SB vine 


Substituting Eqs. (3), (8), and (14) into Eq. (13) above produces 


Elrod's [1] "propagation equation" for stagnation enthalpy: 


1 DH = 7?7H+VH-Vino (15) 
a Oo DE \ 


For spherically symmetric flow in a conical passage with slight 
flow deviations which are also spherically symmetric, the flow is des- 


cribed by only one spatial coordinate - namely the radius r originating 


1] 





at the apex of the cone. 
Chiang, et al., [2] used linear perturbation theory to express a 
quantity as the sum of a steady-state value and a perturbation. 


Stagnation enthalpy was written as 
= ? 


where the prime denotes a perturbation term. 


Since the stagnation enthalpy is uniform throughout the mean flow, 


it is apparent that 


VH=0 (17) 
and 

Dia) = 18 

am: ue 


Hence, Eq. (15) may be written, to the first order in small quan- 
tities, as 
D || DH’ VA q 
ee |) +VWH > Ven (19) 
Dt |Q* Dz C 
As mentioned previously, spherically symmetrical flow is one- 


dimensional. Chiang, et al., [2] pointed out that for this case 


= 9 
or 
2 ECS alee: a 
Va, 2(r 2) (20) 
oe + GO 
Dee OL or 


The upper sign pertains to diverging flow and the lower sign is 
related to converging flow. This choice of signs is based upon the 


fact that for diverging flow the velocity and increasing radius are 
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identical in direction. For converging flow, the directions are opposed 


Therefore, the lower (negative) sign applies to nozzle flows. 


For convenience, Elrod [1] introduced the following dimensionless 
variables: 





—~ H’ 4 — O.,% : = 
p= az ) mw ) t R (21) 


where ay is the acoustic velocity at stagnation conditions and R is 


the fictitious radius at which M = 1 for conical flow. In terms of 


these new variables, the "propagation equation" Eq. (19) becomes 


Se[le) Se] --& 39 [(G) 3 a 


where 


D- 2 + Gd 
DT oY Go a¥ a 


Elrod [1] considered sinusoidal disturbances and found it bene- 


ficial to perform a frequency analysis by representing in the form: 


6=F($)e"™ 


(24) 
W is a dimensionless frequency number. 
It was also assumed that 
F(<)= F-)* exp friwfmas \ n(¥) ~ 
The resultant differential equation for n(&) is 
n’(é) +( W?- #(&) | A(§) =O (26) 
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where 


foe (27) 
35 |& dy (28) 
and 
Ac a \2 mae 
g(t) =~ a nts ahd aol) 
- (29) 





ee 5 E-3. y-) 37-5 
= M 
| -ae>o | 


Using the steady state, one-dimensional, isentropic relationships 
for area and temperature as functions of Mach number (see Shapiro [5]), 


the following relationships hold 


- 
pier i.e. || +a! Me (30) 
Be =(te) = [I st 


wt! 


io ; See 
= =(A,\". “e{ 2 (f+ xt m2) |4O 
Y= R =(4 i E | (31) 


+ Chiang [2] arrived at the general solution for Eq. (18) by using 
the WKB method. (See Morse and Feshback [6]). 
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AE) = Cewpl-iwe do fale) ag | 


+C,exp[awg-dq[a@dg] 
C, and C., are arbitrary constants. £q. (24) now is written as 
g=(* exp| * Mes ¥ cieap ew(T-¥)] 
x exp or J 9) de) + Cz exp iw(r+$)| 


x enp|-dy ( ale) 43) 


(33) 


Chiang [2] states that the C, term with the exp [ iw(t-e)] indi- 
cates an outgoing wave (in the positive r-direction) while the C, term 
with exp [ iw(t+8) | indicates an incoming wave; the term exp [ iw f mde | 
Peace cetits a mean convection; and the factors exp Ea g(e)de | 
represent attenuation of waves. For the case of an infinite conical 
nozzle, with pressure fluctuations impressed at the throat, C, = Q. 

The WKB solution is asymptotically valid in the limit of small Mach 
number M and high frequency number W with an error of the order (MP /w°). 
Chiang [2] extended the solution to larger Mach number and smaller 


a“ 


frequency number W. A summary of this extension follows: 


LS 





Assume 


o Law 


n= 6 EZ (8) +e, FZ 08) on 


After substituting into the previous differential equations for n 


and making comparisons, asymptotic approximations of Z,(E,) and 


nah f wt - (35) 

a a 
"9(8) 4 3 

= 1-4 | ole) de’ 


where 55 is sufficiently en far enough away from the throat) 


Z,(&0) are determined. 


z,(¥) 


THe 


2(9 


[f 
s ae 


(36) 


so that the WKB method holds. Differentiating Eqs. (35) and (36) 
yields 


z/(é,\=s5 H%.) < 


ee, 3.) (38) 


It is further shown that Z, and Lo are complex conjugates and that 
Eqs. (35) - (38) can be used for initial conditions to solve the 


differentia! equations for Z, and Z,. 
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Feet 


After changing variables by using 


X= M* 


and 


Chiang showed that 
$= -4+484 
iy 
and arrived at the following differential equations 
fod +] SE aso + aw w/(] 
_ 4 Po) =O 
V(x) Sake W'(x) -2W 4] 


~ 33) U- (x) =a) 


(39) 


(40) 


(41) 


(42) 


(43) 


(44 ) 


(45) 





with the initial conditions 


U,(X.) =! 

Uy (%) =O i 
Xo 

stein (a 


where Xo is the value of x corresponding to &€ = - Bea Sali eallotalae 
have all been previously defined as functions of M and are therefore 
functions of x. 

It is always possible to select a Mach number M small enough to 


Satisfy the restriction on the WKB solution, 


3 
ae << | (47) 


The amplitude and phase relations between pressure and mass 
flux perturbations are determined in the following manner. Eqs. (8) 


and (10) are used to obtain 


29! - 7" OF | (48 ) 





Ox or 
ies. 8H’ DH” i 49 
wae Ge Cf Or _ 
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The mass flux perturbation is written 


(ea) = et oe (60) 


After introducing the dimensionless pressure and mass flux perturba- 


tions: — 


(51 ) 


Be (04%). | | (52) 


Chiang [2] showed that 











oT dind+M_ 2ing 

Sy Meee (53) 
7 m 2A <= Ln 

3g ay a 


Furthermore, since both m and pu are periodic in time and contain the 


factor exp [ iwt | 








olT 
JT = 27 (54) 
ye a7 


For the case of an infinite conical nozzle with sinusoidal 
pressure fluctuations impressed at the throat, C, in Eq. (34) be- 


A 


comes zero and 


19 





(8) = C exp |-LW é | (u,+2uy) 
b= C, (s* cxpl-hw [Md | 
x exp[iw(r- 9) ] G, +ivy) 


Taking the logarithm of @ and differentiating gives 


oT 
pecoere (i = ¥-1 db Au 
ag oly gt) + cea (4 Ax 


a 4y)\ " Lh —w xe 


+ Ban (4 & - 44) 


Substitution of Eq. (57) and (58) into Eq. (53), together with 
Eq. (54) yields 





ol 
eon, . 1M 7 e** 
fae ra 
yo ee 1-™ 
where 
Cae Nie = 
— + 
NV = amplitude factor = LS 
NS tM, 
| WY ss phase angle = arctan [* N,- N,N, 
: N, No te N. 
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(55) 


(56) 


(57) 


(58 ) 


(59) 


(60) 


~ (61) 








pl- 


N, = —N, -W% 


A dimensionless parameter designated as the modified Strouhal 
number S* was defined by Chiang [2] 


i¥ Gant 
s¥ ao" 


oe (63) 


and shown to be related to the frequency number 
at 
_ MS 
Wee 


For a given Mach number M, it is more convenient to express the 


mass flux-pressure relationships in terms of the modified Strouha!l 


number S* rather than the frequency number W. 
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III. NUMERICAL SOLUTION 


Equations (44) and (45) are transformed to a system of four first 


order differential equations by defining the variables: 


yz (x) = uy (x) 
yy(x) = uy (x) 
2 T (65) 
Y3 (x) = vy (x) 
Yq (x) = vy (x) 
The resultant system of differential equations is 
/ 
yy * Yo. 
pan dP 4 +2Wy4] + GAZ 
¥9 = 2 lax @ 7 | Ee al 
(66) 
YS= Ye 


, (x) 
egg en th] 92 


A numerical solution for this system of equations was obtained 
using Hamming's Modified Predictor-Corrector Method of numerical inte- 
gration. (See Ralston and Wilfe [7}). It is a stable, fourth order 
integration technique which uses four preceding points for the evalu- 
ation of each successive point. PROGRAM TWO makes use of an IBM [8] 
Subroutine, DHPCG, to perform the numerical integration. 

Evaluation of the initial conditions V1 (x) and V1 (x0) is accom- ° 
plished using PROGRAM ONE. Subroutine GOPHI evaluates the function 


g(x)/o. Numerical integration of V1 (x,) is performed using a twelve 
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~—.~ 


point Gaussian Quadrature formula which integrates polynomials up to 
degree 23 exactly. Subroutine DQG12 is used for this purpose. (See 
IBM [8], p. 299.) 

With uw), u,, v,, and v, known, subroutine OUTP of PROGRAM TWO 
computes values of amplitude Factor N and phase angle w for a given 
Mach number M and frequency number W. The corresponding modified 
Strouhal number S* is also computed. 


The programs have been written for the case of air with y = 1.4. 


IV. EXPERIMENTATION 


A. EXPERIMENTAL APPARATUS 

To approximate the characteristics of an infinite conical nozzle, 
a nozzle 17.5 inches in length with an 8 inch diameter entrance and a 
1 inch diameter exit was designed and constructed. It subtended a 
half-cone angle of 11.3 degrees and featured an entrance-exit planar 
area ratio of 64. This ensured that stagnation conditions existed at 
the entrance for all subsonic flows in the nozzle. (See Shaprio [5], 
p. 614.) 

Figures 1 and 2 show the overall arrangement of the experimental 
apparatus and associated electronic equipment. Figure 3 schematically 
depicts the test section. 

Compressed air from two supply tanks passed through a pressure 
regulating valve and into the 10 inch diameter stagnation chamber 
through a baffling arrangement. The stagnation chamber contained 
honeycomb material with a length/diameter ratio of 10 to dampen out 
any disturbances in the plenum air. A fine mesh screen was also installed 


to hasten turbulence decay prior to the entrance of the flow into the nozzle. 
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Sinusoidal pressure disturbances in the nozzle exit plane were 
created by a vibrating piston located ‘nethe exhaust chamber. The 
piston was driven by a MB Electronics Model EA 1250 Vibramatic Exciter 
powered by a Model 2120 Amplifier using a sinusoidal signal froma 
Hewlett-Packard Model 200CD Wide Range Oscillator. The vibrator was 
mounted on three horizontal tracks so that the zero position of the 
piston could be adjusted for optimum operation. 

Dynamic pressure measurements in the exit plane were recorded with 
a Bruel and Kjaer Model 4138 1/8 inch Condenser Microphone with Adaptor 
UA0036 and Model 2618 1/4 inch eae connected to a Model 2803 
Two-Channel Power Supply. 

Mass flux (oq) measurements were performed using a Thermo-Systems, 
Incorporated Model 1050 Constant Temperature Anemometer with a DISA 
Type 55F33 Hot Wire Probe. The hot wire anemometer measured Doth steady 
State (DC) and dynamic (AC) components of mass flux. The anemometer DC 
voltage was measured on a Simpson Model 2700 Digital Voltmeter. 

The anemometer and microphone signals were passed through identical 
Krohn-Hite Model 3750 Filters (band-pass) to preserve phase. The AC 
voltages were measured by two Hewlett-Packard Model 3400A RMS Voltmeters. 
The filtered dynamic signals were displayed on a Hewlett-Packard Model 
132A Dual Beam Oscilloscope for visual verification of system response. 

Phase angles between pressure and mass-flux fluctuations were meas- 
ured on a Dranetz Model 305-PA-3002 Digital Phase Meter using the output 
Signals from the oscilloscope amplifiers. Figure 4 shows a block diagram 
of the electronic measurement equipment. 

Manometers were used to measure stagnation and exit static pressures. 


Stagnation temperature was measured by a mercury-in-glass thermometer. 
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B. PROCEDURE 

Mass flux (pq) calibration of the hot wire anemometer was performed 
using a Thermo-Systems, Incorporated Model 1125 Calibrator. The method 
of least squares was used to fit an eighth order polynomial in the 
quantity (v2 - v4) to the calibration points. (See Computer PROGRAM 
THREE.) Microphone calibration was performed with a Bruel and Kjaer 
Piston Phone: Type 4220. 

Each experimental run consisted of first setting a steady state 
flow rate by maintaining a constant pressure in the stagnation chamber. 
The frequency and amplitude of the pressure disturbances were set on 
the oscillator. Both band-pass filters were set to pass only a narrow 


band at the disturbance frequency with gain settings of 20 db. Data 


was taken over the range of frequencies from 20 Hz to 250 Hz. 


C. DATA REDUCTION 

Reduction of experimental data was accomplished ona Seger com- 
puter through the use of PROGRAM FOUR. 

Steady state mass flux was computed by using the anemometer DC 
voltage in the evaluation of the hot wire polynomial. The magnitude of 
mass flux perturbation was obtained by first correcting the anemometer 
RMS voltage and then adding it to the DC voltage to yield a peak voltage. 
‘The hot wire polynomial was then evaluated for the peak voltage, giving 
a maximum value of mass flux. The steady state value was then subtracted 
from the maximum value to arrive at the perturbation. 

Microphone RMS voltage was similarly corrected to obtain a peak value 
and then multiplied by the microphone sensitivity to yield the magnitude 


of the pressure perturbation. A 6 db. narrow band insertion loss in the 
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band-pass filters was also accounted for in the computations of dynamic 
Quantities. 

Mach number was computed from the one dimensional, isentropic 
pressure relationship. The velocity q used in the computation of 
Modified Strouhal number S* was determined by dividing the mass flux 
by a value of p arrived at by using the perfect gas law and the isen- 


tropic pressure-density relationship. 


V. DISCUSSION OF RESULTS AND CONCLUSIONS .. 


Figures 5 through 14 show a comparison of the experimentally deter- 
mined amplitude factors and phase angles with values predicted by the 
Modified WKB solution for Mach numbers: .084, .117, .201, .278, and .340. 
Both amplitude factor and phase angie are plotted against modified 
Strouhal number S*. The smooth curves are the modified WKB solutions 
and the experimental results are represented by crosses. 

For low frequency sinusoidal pressure disturbances with correspond- 
ing small values of S*, the phase angle is -180 degrees in accordance 
with the modified WKB solution which asymptotically approaches that 
value. An increase in pressure results in a decrease in mass flux and 
physically it is easy to visualize an effective phase angle between the 
rising pressure and the falling mass flux. Mathematically, this effec- 
tive phase angle is simply 180° plus the value of the actual phase angle. 
Thus, for low frequencies there is no effective lag between 7m and u. 

The WKB values of amplitude factor monotonically increase from unity 
for low S* to the asymptotic limit for extremely high S*. This means 
that the magnitude of the mass flux response diminishes as the frequency 


of the pressure disturbances increases. 
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Chiang [2] noted that the quasi-steady relationship was probably 
valid for low values of S* but became increasingly inaccurate for larger 
values of S*. 

The experimental values of amplitude factor agree quite well with 
those predicted for all Mach numbers tested. The experimental values 
for phase angle differed from those predicted but exhibited a similar 
trend. For the low frequency tested (20 Hz) all experimental phase 
angles were close to -180 degrees while the WKB values were slightly 
higher. For higher frequencies, the experimental phase angles became 
larger than the WKB predicted values. This apparent phase shift can 
be attributed to effects of viscosity in the flow. Viscosity serves 
to dampen the flow system with a resultant increase in effective phase 
angle over the Wloaeee eee 1.e., that assumed in the WKB analysis. 

Pnysicai Timitations on the test apparatus prevented investigations 
in the regions above M = .34 and disturbance frequencies of 250 Hz. 
Electronic noise and vibrations penenated by the air compressor prevented 
its being run during experiments. Consequently, the air supply in the 
tanks could not sustain runs over a frequency range at higher Mach num- 
bers. The other problem encountered was that of turbulence. Early in- 
vestigations revealed that turbulence in the nozzle was prevalent at 
frequencies below 200 Hz. Inclusion of the fine mesh screen in the 
Stagnation chamber helped greatly to dissipate turbulence or at least to 
break it up into high frequency turbulence. Power density spectrums of 
the hot wire anemometer signals, obtained by taking the Fourier Transform 
of the autocorrelation function, indicated the presence of turbulence 
with frequencies above 250 Hz. This high frequency turbulence created 


too high a signal to noise level for meaningful data to be obtained 
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above 250 Hz. A PAR Model 101A Correlation Function Computer and a 
Model 102 Fourier Analyzer were used to determine the power density 
spectrums. 

It is concluded that the modified WKB solution, with its ideal 
flow assumptions, does offer a reasonable approximation of the ampli- 
tude factor for pulsating subsonic flow of air ina conical nozzle, 
as long as the flow originates at near-stagnation conditions. Pre- 
diction of phase angle is conservative and could be used as a lower 
bound. 

t is recommended that further investigations be pursued at higher 
Mach numbers. This would require both a quiet compressed air supply 
and a higher power vibrational exciter. It is felt that the present 
location of the air compressor within the laboratory is unsatisfactory 
and chat the compressor snouid be isoiated Trom tne iaboratory by 
either moving it outside of the building or by some other suitable 
means. 

It is also recommended that an analysis be undertaken which would 
account for the effects of viscosity upon the dynamics of a pulsating 


Flow field in a convergent nozzle. 
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Figure 1. TEST APPARATUS 
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Figure 2. ELECTRONIC EQUIPMENT 
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FIGURE 5. COMPARISION OF EXPERIMENTAL AMPLITUDE FACTOR WITH 
MODIFIED WKB SOLUTION, M = .084 
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FIGURE 6. COMPARISON OF EXPERIMENTAL PHASE ANGLE WITH 
MODIFIED WKB SOLUTION, M = .084 


ge 


-90 





30 





10 
1.0 
S* 
as 
0.1 
703 
S555 20 220 S40 4.0 5.0 6.0 


AMPLITUDE FACTOR N 


FIGURE 7. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .117 
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FIGURE 9. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .20] 
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FIGURE 10. COMPARISON OF EXPERIMENTAL PHASE ANGLE 
WITH MODIFIED WKB SOLUTION, M = .201~ 
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FIGURE 11. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .278 
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FIGURE 12. COMPARISION OF EXPERIMENTAL PHASE ANGLE 
| WITH MODIFIED WKB SOLUTION, M = .278 
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FIGURE 13. COMPARISON OF EXPERIMENTAL AMPLITUDE FACTOR 
WITH MODIFIED WKB SOLUTION, M = .340 
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FIGURE 14. COMPARISON OF EXPERIMENTAL PHASE ANGLE ~ 
| WITH MODIFIED WKB SOLUTION, M = .340 
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